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Abstract 

The longitudinal-transverse spin asymmetry Alt in the polarized Drell-Yan pro- 
cess depends on the twist-3 spin-dependent distributions of nucleon. In addition to the 
contributions expressed as matrix element of the twist-3 operators, these distributions 
contain the so-called Wandzura-Wilczek part, which is completely determined by a 
certain integral of the twist-2 spin-dependent parton distributions. We demonstrate 
that the recently obtained empirical information on the transversity distribution al- 
lows a realistic estimate of the Wandzura-Wilczek contribution to Alt for the case 
of polarized proton-antiproton collisions. In particular, our results in the Wandzura- 
Wilczek approximation indicate that rather large Alt can be observed in the proposed 
spin experiments at GSI, and its behavior as a function of dilepton mass obeys novel 
pattern, compared with the other double-spin asymmetries Att and All- Our results 
provide a guide for testing a signal of effects originating from the twist-3 operators 
associated with quark-gluon correlation. 
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The proposed polarization experiments with antiprotons at GSI [T] stimulate renewed 
interest in the polarized Drell-Yan processes to access the chiral-odd spin-dependent parton 
distributions of the nucleon. The double transverse-spin asymmetries A-tt for lepton pair 
production in collisions of transversely polarized protons and antiprotons, p^p^ —>■ l~^l~X, 
are estimated, and are found to be large enough to be measured at GSI [2j, providing promis- 
ing way to probe the chiral-odd twist-2 spin-dependent parton distribution, the transversity 
hi{x) [31 m [5l [6l [7] . In particular, the pp collisions at moderate energy in GSI experiments 
allow us to probe the relevant parton distributions in the "valence region", in contrast to 
the complementary case of pp collisions at, e.g., RHIC where the "sea-quark region" is 
mainly probed. The QCD corrections to Att at GSI have been studied recently at next-to- 
leading order (NLO) [8] and at higher orders with the "threshold resummation" [9]. The 
resummation corrections relevant when the transverse-momentum of the produced lepton 
pair is small [TOl [H] are also investigated [12]. It has been found that the behavior of 
these QCD corrections associated with the valence region is rather different from the cor- 
responding effects involving the sea quarks for the pp-collision cases [131 [HI [H]. As a 
result, these QCD corrections are small at the kinematical regions corresponding to the 
GSI experiments, suggesting that the large LO Att at GSI is rather robust. This fact also 
allows us to estimate the value of Att at GSI using only the empirical information on the 
transversity distributions [12], which is recently extracted [15] through the LO global fit 
to the semi- inclusive deep inelastic scattering (SIDIS) data, in combination with the e'^e~ 
data for the associated (Collins) fragmentation function. 

The PAX Collaboration has proposed the Drell-Yan experiments in pp collisions at 
the CM energy ^/s with s = 30 and 45 GeV^ in the fixed-target mode, and those up to 
s = 210 GeV^ in the collider mode [1]. Those GSI- PAX experiments will measure Att 
for 0.2 < Q/y/s < 0.7 with Q the mass of the produced dilepton, and indeed probe the 
transversity hi{x) in the valence region in a wide range of x. It should not be overlooked 
that the double-spin longitudinal-transverse asymmetry Alt is also readily accessible in 
those Drell-Yan experiments, in particular, in the fixed-target mode with the longitudinal 
polarization of the target: Alt plays a distinguished role in spin physics because it allows 
us to access the twist-3 spin-dependent parton distributions as leading effects [1], similarly 
as the longitudinal-transverse asymmetry associated with the structure function g2 in the 
polarized DIS [16]. Thus the data of Alt will provide an experimental test whether the 
quark-gluon-quark correlations inside the nucleon is sizeable or not, in particular, in the 
chiral-odd spin structure that is not accessible by g2 in DI^. These facts call for theoretical 
study of Alt to assess its potential at GSI experiments, which is the purpose of this Letter. 
Up to now Alt was estimated for the pp-collision cases [18], but the above mentioned 
situation for Att suggests that the behavior of Alt will be different between the pp and pp 
collisions. Also it is important to clarify the impact of the new empirical information [13 [IS] 
of the transversity distribution hi{x) on the prediction of Alt at GSI, because Alt depends 
on hi{x). We will demonstrate that this empirical information for hi{x) indeed allows a 
useful estimate for Alt in PP collisions at GSI kinematics. 

*For a test of chiral-odd quark-gluon-quark correlation using the SIDIS data, see [T7]. 
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We shall work at LO QCD, which provides a sufficient accuracy for our first estimate 
of Alt at GSI. We may anticipate that the mechanism associated with the valence region 
relevant to GSI kinematics could make the QCD corrections to Alt small, similarly to the 
case for Atx mentioned above. To calculate A^t, "we first recall the parton distributions of 
the nucleon. At LO, we need the spin-dependent quark distribution functions of twist-3 as 
well as of twist-2, which are defined as the nucleon matrix element of the chiral-odd and 
chiral-even quark bilocal operator with the light-like separation between the constituent 
fields [lEllTl, 

= 2 [h{x, /i^) {S±,,P, - S^,P,) /M + hL{x, /i2)M {P^n, - P,n^) {S ■ n)] , (1) 
/ ^e^^^PS\mi,l6i^i>^n)\PS) = 2 [g,ix, fi')P,{S ■ n) + g^x, fi')S^,] , (2) 

and we also need the unpolarized quark distribution defined as usually as 

I ^e''^{PS\miM^^)\PS) = 2/i(x,/i^)P^, (3) 

where \PS) denotes the nucleon state with mass M, the four momentum P^ = {P^, P~ , 0±), 
and the spin vector 5^ satisfying P^ = 2P+p- = M^, 5^ = -M^, and P ■ S* = 0, and a 
light-like vector = (0, n~, 0^) is introduced by the relation P-n = 1. is decomposed as 
S^' = {S ■n)P^-M^{S -71)71^ +S'[ with P-S± = n-Si_ = 0. In (P)-©, the gauge link operators 
which ensure gauge invariance are suppressed for simplicity. The distribution functions /ii,L, 
gi^T and /i depend on the factorization scale /i, at which the bilocal operators in the LHS 
are renormalized, and those distribution functions are defined for each quark and anti-quark 
flavor ip = ip^ [a = u, u, d, d, s,s, . . .) as etc. We remind that in the infinite momentum 
frame {P^ oo) the Lorentz structures associated with hi, gi and /i are of O(P^) (twist- 
2), those for hi and gx are of 0(1) (twist-3), and those behaving as twist-4 (0(1/P+)) 
are ignored in the RHS of ([I])-®. g\{x,^^) is the familiar helicity distribution for u-quark 
carrying the momentum component = xP^ inside the longitudinally polarized nucleon, 
and, similarly, /i"(x, /i^) is the transversity distribution inside the transversely polarized 
nucleon [3l H]. Note that the twist-3 distributions hi and gr are also associated with the 
longitudinal and transverse polarization of the nucleon, respectively. 

The above mentioned classification of twist based on the power counting in the infinite 
momentum frame is directly related to the power of 1/Q with which the corresponding 
distributions appear in the physical cross sections, but does not exactly match the conven- 
tional and formal definition of twist as "dimension minus spin" associated with the relevant 
operator structure in ([T]) and ([2]). As a result, the distributions and gx actually contain 
the piece that is expressed by matrix element of the twist-2 operators as [201 IH El (see also 



3 



Appendix in [21] ) 



hlix^^i') = 2x rrfy^%^ + ---, (4) 

J X y 



gUx,f^') = I + (5) 



1 ^at„^ ,,2\ 

y 

where the eUipses stand for "genuine twist-3" contributions given as matrix element of the 
twist-3 operators; it is known that those twist-3 operators can be reexpressed as quark- 
gluon-quark three-body correlation operators on the lightcone, using the QCD equations of 
motion [H [221 [23l [2l] . In the following we call the twist-2 component, shown explicitly in (j4]) 
and ([5]), the Wandzura-Wilczek part. Because the operators with different geometric twist 
do not mix with each other under renormalization, the Wandzura-Wilczek part does not mix 
with the genuine twist-3 contributions under the QCD evolution with /j?. Thus both x- and 
//^-dependences of the Wandzura-Wilczek part are determined solely by those of the twist-2 
distribution functions as (j4]) and ([5]). Taking into account only the Wandzura-Wilczek part 
in (jl]) and ([5]) yields the 'Wandzura-Wilczek approximation" for and qt- 

With the above definitions for the parton distributions, we can write down the LO 
expression for the longitudinal-transverse spin asymmetry Alt in the polarized pp collisions. 
Before doing this, it is worthwhile to remind the LO formula of the other double-spin 
asymmetries A^l and Att OIHE]. Using the quark distributions inside the proton, 

< dQ^dxpdfl dQ^dxpdfl 



dcr^^ 






dQ^dxpdQ 




dQ'^dxpdQ 


da^^ 


+ 


dcr~^'~ 


dQ'^dxpdQ 


dQ^dxpdU 








dQ'^dxpdQ 




dQ'^dxpdQ 


dcrTT 


+ 


daU 


dQ'^dxpdO, 


dQ^dxpdft 



A at^-axpaiL at^'-axpaa ^ ^a^l 

■ ''^''Y.a^in{x,,Q^)n{x,,Q^y ^ > 

for the production of the dilepton with the invariant mass Q and the longitudinal momentum 
component Qz corresponding to the Feynman xp-, where one of the leptons outgoes to 
the direction with the angle VL = {6,(f)). Ca represents the electric charge of the quark- 
flavor a and the summation is over all quark and anti-quark flavors, a = u,u,d,d, s,s, . . .. 
The scaling variables xi^2 represent the momentum fractions associated with the partons 
annihilating via the Drell-Yan mechanism, such that = [xiPi + X2P2Y = X1X2S and 
xp = Xi — X2 {= 2Qz/ y/s in the CM frame), where s = {Pi + P2)^ is the CM energy squared 
of the colliding proton and antiproton. This implies 



- . / 2 4Q2 \ 1 / / 2 

= - \xF + yxp + I , ^2 = 2! + yxp + I . (8) 

In (|6]) and ([7]), d^L and d^r represent the asymmetries in the parton level deflned as 

sin^6' cos2</) , . 

a-LL = 1, orr = -r- 5^, (9) 

1 + cos'^o' 
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with the polar and azimuthal angles 6 and (p in the dilepton rest frame with respect to 
the incoming beam and transverse-spin axes, respectively. The LO formula for A^t in 
p~^p^ X or p^p~^ l^l~X can be expressed similarly as [U [6] 



dQ'^dxpdQ 


dQ'^dxpdfl 


da^^ 

.1^0 7 . ir-\ 


1 da^^ 



A dQ^dXpdM dQ-^dXpdil ^ 

^LT = T— ^ T— T = aLT 



2^ 



(10) 

associating the variables xi and X2 with the longitudinally and transversely polarized beams, 
respectively, with 

M 2 sin26' coscj) , , 

= Q 1 + cos^^ • 

We note that A^^l and A'yx receive contribution only from the twist-2 distributions, while 
Ai^T is proportional to the twist-3 distributions and hence ai,T is suppressed by a factor 
l/Q compared with ([9]). 

To compute the above formulae ([6]), ([7]), and (fTOj) with the GSI kinematics, we have 
to specify the LO parton distributions to be substituted. We use the LO GRV98 [25] and 
GRSV2000 ("standard scenario") [26] distributions for the unpolarized and longitudinally- 
polarized quark distributions /f (a;, Q^) and (x, Q^), respectively. For the LO transversity 
distribution h\[x^ Q^), we are guided by the recent information from the LO global fit [151 
[T9] : we find that a useful estimate can be obtained by assuming the relation 

h\{x,^^)=gl{x,^^), (12) 

at a low scale \i (/i^ = 0.26 GeV^ using the GRSV2000 gi{x, fi"^)); its QCD evolution from 
IJ? to is controlled by the LO DGLAP kernel [27] for the transversity. It is worth noting 
that the above relation (fT2|) at the low /i^, which is exact in the non-relativistic limit, is 
suggested also by the estimates from relativistic quark models for nucleon [6], [2], [28] , matches 
the results by lattice QCD simulation [29], [30], and has been used in the previous estimates 
for Att at GSI [21 [HI [12]. The obtained LO transversity distributions for u and d quarks, 
xh^{x,Q'^) and xhf{x,Q'^), are shown in Fig. [11(a) as a function of x a.t Q"^ = 2.4 GeV^. 
xg^''^{x, Q"^) is also shown in the same figure. For convenience, we have multiplied the factor 
— 1 to the (i-quark distributions. If we compare Fig. [H^a) with the results of the LO global 
fit [miin], we see that, for the valence region 0.2 < x < 0.7 relevant for the GSI kinematics, 
our LO transversities lie slightly outside the error band of the fit, similarly as observed for 
the NLO case [12]. Therefore, our transversities will provide a realistic estimate of the 
upper bound of the relevant asymmetries, implied by the present empirical uncertainty in 
the transversities. (At present there are no data to constrain the transversity hi{x,Q^) 
directly for x > 0.4, and in this region the uncertainty bands resulting from the LO global 
fit [151 [IS] could be subject to the particular choice of the parameterization of h1{x,Q'^) 
assumed in the fitting procedure.) We see from Fig. [TJ^a) that (/i"(x, Q^))^ ^ {hf{x, Q^)) 

and {gi{x,Q'^))'^ ^ {gi{x,Q'^)Y in the valence region, and likewise [23] for /"(x,(5^). This 
also holds for higher Q^, so that 

Att _ hUxi,Q')hUx2,Q') . . 

aTT mx^,Q^)mx2.Q^y ^ ^ 
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Figure 1: (a) The transversity distributions xhi{x,Q^) and the hehcity distributions 
xgf{x,Q^) at the scale = 2.4 GeV^ for u- and d-quarks. (b) The twist-3 distributions 
xh\{x,Q'^) and xg^{x,Q'^) in the Wandzura-Wilczek approximation for u- and (i-quarks. 
For convenience, we multiphed —1 for d-quark distributions in both figures. 

for ([7]) at GSI, and hkewise for Hence the GSI measurement of ([7]) allows a direct 
access to hl{x,Q'^) P 181 fT2f| . 

The twist-2 spin- dependent distributions obtained above determine hi and qt in the 
Wandzura-Wilczek approximation using (jl]) and The results for the u- and d-quarks 
with = 2.4 GeV^ are shown in Fig. [T](b), with the curves for the d quark showing the 
results multiplied by —1. The integral with the factors 1/y, in (jlj), ([5]) shifts the 
peak of the curves to lower x with the suppressed peak-height, compared with those for the 
corresponding twist-2 distributions in Fig. [D^a). We also see the w-quark dominance in the 
valence region, similarly as in the twist-2 distributions, so that 

for ffTOj) at GSI kinematics, where the ellipses denote the contributions associated with the 
genuine twist-3 operators. This implies that the 'Wandzura-Wilczek contribution" to Alt 
at GSI is directly related to the behavior of /i" (x, Q^). 

t For the values of Q > 1 GeV relevant to the Drell-Yan process, our LO transversities satisfy Soffer's 
inequality [31], 2\h'{{x, < fi{^i Q"^) + dii^^ Q^)^ for the u- and d-quarks, but violate it for the d-, s- 
and s-quarks by a small amount, similarly as in the previous works [12|ll8j. because of negative polarization 
for the d-quark, gi{x,Q^) < 0, and the smallness of s- and s-densities (/f(a;,/i^) = /f(a;,^^) = at the 
input scale n in GRV98). This violation of Soffer's inequality will be harmless to our numerical estimates 
of the asymmetries because of the dominance of the u-quark distribution as (jl3[) and (jl4p . 
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Figure 2: (a) Alt and (b) All and Att as a function of xp for Q = 2.5 and 4 GeV at 
s = 45 GeV^. For A^t-, the chiral-even and -odd contributions are also shown separately. 



The data on the transverse spin structure function g2 from the polarized DIS experiments 
indicate that the genuine twist-3 contribution in ([5]) is small [I6] , and qt approximately 
follows the Wandzura-Wilczek result. The calculations of low moments of qt by lattice QCD 
simulation support this result [301132] . Also, estimates from nucleon models, combined with 
the QCD evolution for the relevant twist-3 operators [331 121] , suggest that the Wandzura- 
Wilczek part of ^ and ^ dominates qt and /il for /i^ > 1 GeV^ [Ml [35l [18] . For the 
present first estimate of Alt of (ITOl) in pp collisions, we employ the Wandzura-Wilczek 
approximation of Fig. (H^b) for qt and Hl- 

In all the following numerical evaluation, we present the results for the "reduced asym- 
metries" Ayw = Ayw/avw W = L,T)§. We first consider the fixed-target mode, where 
Alt will be readily accessible. Figure [2t^a) shows Alt as a function of a;^ for Q = 2.5 and 4 
GeV at s = 45 GeV^. Also shown are the separated contributions from the chiral-even and 
-odd distributions, corresponding to the first and second terms in the numerator in ffTOj) . 
The results may be compared with the behavior of All and Att at the same kinematics, 
shown in Fig. [2](b). The curves for Att reproduce the corresponding LO results in [S]. All 
and Att are symmetric with respect to xp = 0, while Alt is not symmetric (compare 
([7]) with ffTOl) ). All and Att are almost flat as a function oi xp for the GSI kinematics [8], 
in strong contrast to Alt- These features of Alt come from the xp dependence of chiral- 
even and -odd contributions; in particular, the chiral-odd contribution shows the tendency 
to increase for decreasing xp, while the chiral-even one shows opposite tendency. The val- 
ues of All and Att are more than 30% and are much larger than their typical values in 
the pp-coUision cases [ISl HE] . This is because of the fact that for the GSI kinematics the 



■Note that qlt of pT|) is defined absorbing the suppression factor M/Q specific to twist-3 cross section. 
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Figure 3: (a) Alt and (b) All and Att as a function of xp for Q = 2.5, 4, 6 and 8 GeV at 
s = 210 GeV^. For Alt with Q = 2.5 and 6 GeV, the chiral-even and -odd contributions 
are also shown separately. 

valence contributions are dominant both in the numerator and the denominator of (El) and 
(ITjIl and the small-x rise of sea-distributions is absent in the denominator [21 [9l El [IS]. We 
see in Fig. |2]^a) that the similar mechanism leads to the significant value (> 10%) also for 
Alt- In general, Alt is smaller than All, ^tt by the presence of the additional factor, xi 
or X2, in (fTOj) compared with (I7j). Further suppression effect for Alt could be caused 
by the behavior of qt and Hl observed in Fig. [U^b) in comparison with Fig. Hj^a). When 
the sea-quark region is probed in pp collisions, these effects, in particular the additional 
Xi^2 factor, lead to Alt much smaller than the corresponding All, ^tt, as demonstrated 
in'[l8]. 

Actually, the fixed-target mode discussed above mainly probes the region xi^2 ^ 0.4 
(see dHl)), where the transversities involved in (fT3|) . (fT^ are poorly determined at present 
(see the discussion above f|T3l) ). In the collider mode we probe the smaller xi^2'- Fig- 
ure [31 is same as Fig. [21 but for Q = 2.5, 4, 6 and 8 GeV and s = 210 GeV^. We 
observe the similar pattern as in Fig. [21 except that in Fig. each of chiral-even and 
-odd contributions changes its behavior between Q = 2.5 GeV and Q = 6 GeV. Also, 
all the asymmetries become somewhat smaller for higher energy, i.e., for smaller Q/y/s. 
Actually, the mechanism relevant to this latter point leads to the behavior commonly 
observed in Figs. [2l and [3l i.e., the increasing All and Att for increasing Q, and the 
corresponding moderate increase of Alt- The corresponding behavior is also presented 
in Fig. m where the relevant asymmetries with xi = X2 = Ql\fs [xp = 0) are plotted 
as functions of Q for s = 30, 45 and 210 GeV^. As clarified in [12], the Q dependence 

^All is slightly larger than Att, because is slightly larger than ft." as in Fig. |T{a). 
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Figure 4: (a) Alt at = as a function of Q for s = 30, 45 and 210 GeV^. (b) All and 
AxT at a; i? = as a function of Q for s = 30, 45 and 210 GeV^. 

of All and Att shown in Fig. ID^b) directly reflects the x dependence of the correspond- 
ing distributions: (fT3l) imphes that Att is controlled by the ratio h^{x,Q'^)/ fi{x,Q'^). It 
is straightforward to see that the scale dependence of the w-quark distributions in this 
ratio almost cancels between the numerator and denominator in the valence region rel- 
evant at GSI, as h'l{x,Q'^)/ f^{x,Q^) ~ /^^'(x, lGeV2)//i"(a:, IGeV^) (see Fig. 3 in [12]). 
Thus the behavior of IGeV^)//^^, IGeV^) function of x directly determines 
the Q-dependence of Att with x = Q/y/s. The same logic holds for All- In the present 
case using GRV and GRSV parameterizations, the ratio hi{x,lGeV'^)/ fi{x,lGeV'^), as 
well as gi{x, lGeV2)//{'(a;, IGeV^), is actually an increasing function of x, leading to the 
Q-dependence in Fig. Ht^b). Note, this mechanism characteristic for the GSI kinematics 
survives even when including the higher order QCD corrections [12]. For Alt, however, 
the cancellation of the scale dependence between the numerator and denominator in (1141) 
is less complete due to the additional y-integral for the Wandzura-Wilczek part, which, 
combined with the additional factor xi or X2 (= Q/y/s), results in the novel Q-dependence 
in Fig. 111(a). In particular, the suppression in the moderate x-region observed in Fig. [T](b) 
compared with Fig. [Tj^a) leads to the decreasing behavior of Alt for increasing Q in the 
large Q region, while the increasing behavior of Alt in the small Q region is caused by that 
of the additional factor 0:1,2 = Q/\fs- 

To summarize, we have presented a first estimate of the longitudinal-transverse spin 
asymmetry Alt for the polarized Drell-Yan process in pp collisions at GSI kinematics. 
Guided by the new empirical information of the transversity, we performed the LO calcula- 
tion of the Wandzura-Wilczek contribution to Alt-, which is directly related to the behavior 
of the transversity in the valence region. The results turned out to be significantly large. 
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and exhibited distinguished behaviors compared with the twist-2 asymmetries Att and 
All. These results serve as a useful guide for possible future A^t measurement at GSI. 

In relation to the discussion on A^x, we also emphasized that the large value of Att at 
GSI kinematics is known to be quite stable when including the QCD corrections, and that 
the behavior of Att as a function of dilepton mass is controlled by the x-dependence of 
the transversity. Thus, first of all, the measurements of Att at GSI will provide the data 
that constrain the detailed shape of the transversity in the valence region, including the 
large x regime where our knowledge on transversity is poor at present. The corresponding 
new information on the transversity will enable us to update our prediction of Alt in the 
Wandzura-Wilczek approximation. If the strong deviation from our updated results were 
observed in the GSI measurements of Alt, this would provide an indication of large genuine 
twist-3 effect, associated with the chiral-odd distribution Hl- The QCD analysis of such 
data using the evolution equation for the corresponding twist-3 operators will reveal the 
quark-gluon-quark correlation inside the nucleon. One problem for this purpose is that the 
exact form of the evolution equation governing the genuine twist-3 contributions in /i^, is 
known to be quite sophisticated even at the LO level [23]. Fortunately, as in the case for 
the similar problem in the chiral-even distribution qt [33], it is proved [21] that, in the 
limit of large number of colors, Nc —* oo, the corresponding evolution equation is simplified 
into the evolution of usual DGLAP-type, with the novel anomalous dimension known in 
analytic form. Since this simplification holds up to the corrections of 0{1/N^) ~ 10%, 
the large- iVc evolution for the genuine twist-3 contributions in Hl provides a powerful and 
practical framework to solve the above problem. 
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